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\  ;  Abstract 

1ft  is  well-known  that  the  complexity  of  performing  operations  on  a  set  depends  heavily  on  the  structure  which 
we  are  allowed  to  put  into  its  representation.  For  example,  searching  through  a  sequence  of  numbers  can  be 
performed  more  efficiently  if  the  numbers  appear  in  sorted  order.  In  this  paper,  we  take,  as  a  case-study,  the 
class  of  problems  involving  a  simple  N-gon  P  and.  making  the  assumption  that  in  addition  to  the  usual 
description  of  the  boundary  of  P,  an  arbitrary  triangulation  is  also  available,  we  investigate  the  computational 
power  gained  from  having  this  additional  information.  Among  other  results,  we  give  a  very  simple,  optimal 
algorithm  for  computing  the  area  visible  from  an  arbitrary  point  hi  P.  We  also  present  several  optimal 
algorithms  for  computing  the  internal  distance  between  two  points  in  P.  Recall  that  the  internal  dhmnca 
between  A  and  B  is  defined  as  the  length  of  the  shortest  path  inside  P  between  A  and  B. 
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1.  Introduction 

The  complexity  of  problems  due  operate  on  fixed  objects  is  highly  dependent  on  the  amount  of 
preprocessing  allowed  in  the  objects’  representation.  As  illustrated  in  the  well  known  paradigm 
searching,  vx  sorting,  the  mere  availability  of  an  order  among  keys  cuts  down  the  complexity  of  searching 
from  linear  to  logarithmic.  In  numerical  analysis,  preconditioning  a  sparse  matrix  is  standard  procedure  in 
order  to  facilitate  the  computation  of  its  inverse.  In  general,  the  crucial  issue  is  to  balance  costs  and  gains  of 
preprocessing  so  as  to  optimize  the  overall  performance.  Few  areas  of  computer  science  are  fine  of  this  type  of 
trade-offs  and.  in  particular,  this  concern  is  recurrent  in  computational  geometry,  operations  research,  and  in 
the  study  of  data  structures  or  data  bases. 

The  first  area  mentioned,  computational  geometry,  provides  a  good  example  of  a  structure.  Le.  the  Vonmoi 
diagram,  easy  to  construct  efficiently,  addle  one  of  dm  moat  powerful  tools  at  our  disposal  for  solving  doeest 
point  problems  [5H77J.  Unrelated,  yet  equally  effective  results  have  shown  that  convex  figures  lead 
themselves  to  speedier  algorithms  than  arbitrarily-shaped  objects  {CH80.CD8QX>K81].  Consequently,  an 
attractive  approach  to  handle  non-con  vex  figures  is  to  decompose  them  into  their  convex  pans,  then  apply  to 
these  the  efficient  methods  known  for  convex  objects  ^H80i:S81.GJ78^C78^V80,TO80].  We  pursue  this 
enoesvor  in  oui  paper*  ana  tnvcsugMG  tne  cmgcncc  or  emcient  aigontnins  ror  various  proot&ns,  assuming 
that  in  addition  to  the  usual  boundary  description  of  a  polygon,  an  arbitrary  trianguiation  is  also  available.  It 
is  standard  to  define  a  trianguiation  of  a  polygon  as  a  convex  decomposition  which  does  not  introduce  new 
vertices1  fCJ78|.  For  our  purposes,  however,  we  can  relax  this  definition  and  allow  the  vertices  of  the 
trianguiation  to  lie  anywhere  on  the  boundary  of  the  polygon.  The  only  provision  to  make  is  that  the  toad 
number  of  points  used  in  the  trianguiation  is  linear  in  the  number  of  vertices  of  the  polygon.  Note  that  this  is 
always  true  with  the  standard,  more  restrictive  definition.  We  also  observe  that  it  is  easy,  given  a  convex 
decomposition  of  a  polygon,  to  derive  a  trianguiation  in  linear  time.  It  is  then  apparent  that  it  is  only  for 
simplicity  that  we  choose  to  be  supplied  with  a  trianguiation  rather  than  a  mote  general  convex 
decomposition  of  the  polygon. 

With  this  additional  information  in  hand,  we  are  able  to  describe  a  very  simple,  yet  optimal  algorithm  for 
computing  the  area  visible  from  any  point  inside  a  polygon.  We  also  present  several  optimal  algorithms  for 
computing  tite  internal  distance  between  two  points  inside  a  polygon.  Recall  that  the  internal  distance  is 
defined  as  the  shortest  distance  a  person  might  travel  from  one  point  to  the  other,  while  remaining  within  the 
boundary  of  the  polygon. 


1  i.c..  when  all  the  vertices  in  Uw  dooomponikm  arc  vertices  of  the  polygon. 
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Next  we  introduce  our  notation,  before  proceedinf  with  the  description  of  the  algorithms.  Let  P  be  a 
simple  polygon  with  vertices  in  clockwise  order.  We  assume  the  existence  of  a  triangulation  TatP, 
defined  as  a  set  of  non-overlapping  triangles  whose  union  is  exactly  P,  and  whose  summits  are  taken  in  the  act 
Tile  edges  of  the  triangulation  which  are  not  edges  of?  are  called  interior  edgpL  As  mentioned 
above,  we  may  choose  to  allow  the  summits  to  fie  anywhere  on  the  boundary  of  P,  provided  that  the  total 
number  of  vertices  in  the  triangulation  does  not  exceed  the  number  of  vertices  of  P,  up  to  within  a  constant 
foctor.  In  this  case,  we  may,  for  simplicity ,  rename  the  vertices  of  Pso  that  the  fist  gives  a  clockwise 

description  of  off  the  vertices  appearing  on  the  boundary  of  the  triangulation. 

Observing  that  a  triangulation  forms  the  embedding  of  a  planar  graph,  we  chooee  a  DCEL  representation 
of  this  graph  as  our  basic  working  structure  [MP78J.  Recall  that  a  DCEL  is  simply  a  handy  data  structure, 
obtainable  in  linear  time  from  any  standard  atgacency  representation,  which  in  particular,  allows  one  to 
traverse  the  boundary  of  each  fine  in  clockwise  order  and  fist  die  tees  encountered  on  die  left-hand  side 
during  the  traversal.  Roughly,  to  each  edge  e  of  the  graph  is  assigned  a  6-fidd  node  containing  the  names  of 
the  two  endpoints  in  some  specified  order,  as  well  at  the  two  adjacent  faces  and  the  names  of  each  of  the 
edges  first  encountered  in  traversing  these  tees  in  clockwise  order,  starting  at  the  endpoints  of  e. 

Note  that  several  algorithms  are  available  in  the  literature  for  computing  an  arbitrary  triangulation  of  an 
N-gon.  The  best  performance  achieved  so  te  is  (XMog  N)  time  [GJ78.CH821  but  is  yet  unknown  so  be 
optimal  or  not. 

2.  Visibility  problems 

A  problem  which  arises  frequently  in  graphics  concerns  the  elimination  of  hidden  lines  from  a  two-  or 
three-dimensional  scene  [NS79J.  In  two  dimensions,  the  problem  reduces  to  computing  the  sets  of  points  that 
are  visible  from  a  given  point  inride  a  polygon  P.  Linear  algorithms  for  this  problem  already  exist 
[CH80.EA81].  but  they  involve  complicated  stack  manipulations  which  become  unnecessary,  once  a 
triangulation  is  made  available.  The  problem  can  be  formulated  as  follows: 

Given  a  simple  polygon  Panda  point  M  inside  P.  the  locus  of  points  V  such  that  the  segment  MV 
lies  entirely  in  P  is  a  simple  polygon  V(M).  Compute  a  clockwise  description  of  the  boundary  of 
V(M)(flg5). 

We  can  regard  the  triangles  of  Tas  forming  the  nodes  of  a  graph  <7.  whose  edges  join  the  pairs  of  triangles 
with  a  common  edge  (i.e*  an  interior  edge)  (fig.1).  As  shown  in  Lemma  1.  the  absence  of  interior  tees 
ensures  that  the  graph  <7  is  actually  a  tree. 


polygon  ■  arid io be Ompfeitt only  adjacent  edges  Memo. 
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[  FIGURE  1) 

Flgerel:  The  triangulatkm  T  and  the  dual  graph  G. 

Lemma  1:  (/Rant. 

Proefc  It  suffices  to  show  that  for  any  pair  of  triangles  tvL  in  the  triangulation,  there  exists  a 
unique  path  between  and  /2  in  G.  The  triangle  ix  partitions  P  into  4  parts.  One  is  the  triangle  ^ 
itsdf.  the  others  being  polygons  adjacent  to  the  edges  of  ^  (note  that  some  of  these  polygons  may 
be  reduced  to  a  single  edge).  At  any  rate,  exactly  one  of  the  three  polygons  contains  the  triangle 
Call  U  this  polygon,  letting  a  denote  its  edge  adjacent  to  ^  and  /  be  (he  triangle  of  T  adjacent  to  u 
and  lying  in  U  (figJX  Since  the  triangulatioa  of  P  also  provides  a  triangulation  of  U,  and  its 
associated  graph  </u  is  a  subgraph  of  G,  we  can  see  that  if  there  is  a  unique  path  in  Gu  from ;  to  ty 
there  is  also  a  unique  path  in  G  from  ^  to  Therefore  we  can  prove  the  lemma  by  induction  on 
the  number  of  vertices.  □ 


[FIGURE  2) 

Fipre  2:  Frevfog  i/hof  (/ h  a  tree 


Let  #be  any  huerior  edge  of  foe  triangulatioa.  and  let  hi  be  any  point  inside  F.  Letting  t  denote  the  triangle 
of  T  which  contains  M,  we  cm  define  GfAf.e)  as  die  unique  subtree  of  G  emanating  than  e,  which  does  not 
contain  /  (figJX 

[FIGURE  3| 

P%mu3:  The  subtree  G(M,$ 


We  are  now  in  a  position  to  give  an  algorithm  for  computing  the  visfoitity  polygon  V(Af).  To  facilitate  our 
task,  we  introduce  the  (Unction  VISIB.  defined  as  fofiows*  let  /  be  any  segment  lying  entirely  on  the  edge  c. 
Remove  from  Tall  the  triangles  winch  do  not  belong  to  G and  call  (2  the  resulting  polygon.  We  define 
VISIB  (M,e*)  as  the  pari  of  Q  which  is  vWMe  from  M  through  the  window  e.  More  precisely.  VISIB  (M,e) 
is  the  set  of  pototiM  in  Q  such  that  foe  only  insersectionofAfu  with  the  boundary  of  g  takes  place  ate*  (fig.4). 
Let  ajbx  be  foe  vertices  of  iw  triangle  in  Q  adjacent  to  *  with  e»ab  and  e  »aV.  Wc  define  A  (resp.  B)  as 
cm  menocckjA  or  roc  po(y|Dflii  mm  (im;  whd  bm  miuum  «tc  pmi|  umHijn  mol  (resp.  mo  /•  it  b  now 


YISfB  (Mj*) 


9  e  He*  on  die  boundary  of  P 

then 

ratm  ({«"}) 

Determine  the  points  ejUL 
9  e  lies  between  A  and  B 


V «-  V1SIB  (M.Bc) 

V*-  V  u  VISIB  (M,cA) 

the 

V «-  VISIB  (MyAB) 
return  ( V)  -  (fig.4) 


[FIGURE  4] 

Flgwre  4:  Computing  the  ana  VISIB  (M.e*)L 

To  complete  the  computation  of  V(AO.  it  suffices  to  determine  the  triangle  of  T  where  M  Hes  -  which  can 
be  done  in  Q(N)  time  -  then  apply  the  previous  procedure  with  respect  to  its  three  edges. 

VISIBIIJTY  ( A40 

Let  eye2,e,  be  the  edges  in  clockwise  order  of  the 
triangle  or  T  which  contains  M.  Initially  V(M)  a  0. 
far  i  a  1,24 

begin 

V(W) V(Af)  u  VISIB 


See  an  illustration  in  figJ.  Note  that,  as  described,  the  procedure  reports  die  boundary  of  V(M)  in 
clockwise  order,  except  for  the  rayedges  of  V(Af),  Lc„  the  segments  coll incar  with  M.  which  are  omitted.  A 
single  pass  through  the  list  V(A/),  however,  will  be  sufficient  to  add  the  missing  segments,  and  we  need  not 
elaborate.  Using  a  DCEL  representation  of  the  triangulation  ensures  that  each  recursive  step  can  be  executed 
in  constant  time,  from  which  we  can  conclude: 

Theorem  2 :  Given  a  simple  N- gon  P  along  with  an  arbitrary  triangulation  of  P,  it  is  possible  to 
compute  the  visibility  polygon  from  any  point  M  inside  P,  in  O(A0  time. 


The  main  advantage  of  this  algorithm  is  that  it  avoids  the  complicated  stack  manipulations  of  [CH80]  and 
[KA81J.  The  reader  may  convince  himsclf/hcrsclf/itsclf  that  the  algorithm  could  be  rewritten  without  greater 
difficulty  in  order  to  deal  directly  with  a  more  general  convex  decomposition  (i.c..  without  first  converting  it 


into  a  triaaguladon).  This  may  be  an  interesting  alternative  if  one  is  willing  to  exploit  the  feet  that  searching 
among  the  edges  of  a  convex  polygon  can  be  done  in  logarithmic  time,  using  a  Fibonacci  search-based 
strategy  {CH80.CD80J.  We  would  not  recommend  this  approach  in  practice,  however,  unless  the  size  of  the 
problem  was  particularly  gigantic.  Once  again,  we  leave  substantiating  these  digressions  to  the  attention  of 
the  reader. 

[FIGURES] 

Figure  S:  The  visibility  polygon  V(MX 

3.  Applications  to  internal  diatanca  problems 

3.1 .  Ttw  ear- racing  problem 

What  is  die  shortest  trajectory  of  a  racing  car  on  a  given  circuit?  More  precisely,  the  problem  which  we 
address  in  this  section  can  be  expressed  as  follows: 

Given  a  simple  polygon  P  and  two  arbitrary  points  A  and  B  in  P,  find  the  shortest  path  inside  P 
between  A  and  B(ftg6X 

[FIGURE  6] 

Figured:  The  internal  path  between  A  and  B.  IP(A.BX 

This  shortest  path  is  called  the  internal  pah  between  A  and  B.  denoted  IF(44)i  and  hs  length,  |IF(4E)l  is 
called  the  internal  distance  between  A  and  B  (flg.6).  To  have  a  visual  representation  of  one  can 

imagine  a  rubber  band  inside  P  tightly  stretched  between  A  and  B.  In  [SM77],  Shames  suggests  an  O (/V*) 
algorithm  for  computing  IF (A,BX  The  method  consists  essentially  of  computing  all  pairs  of  vertices  visible 
from  each  other,  in  OUV2)  time,  so  as  to  form  die  so-called  viewability  graph  of  P.  We  next  add  weights  to  the 
graph  by  associating  to  each  edge  the  Euclidean  distance  between  Its  endpoints.  Computing  an  internal  path  * 
is  now  equivalent  to  finding  the  shortest  path  between  two  vertices  of  a  graph  with  N  vertices,  which  can  be 
done  in  O ( W2)  time.  Of  course,  we  assume  in  this  case  that  both  A  and  B  are  vertices  of  P.  Wc  will  next  show 
how  the  use  of  a  triangulation  permits  us  to  compute  the  internal  path  in  0 (/V)  time,  without  even  having  to 
restrict  the  points  to  be  on  the  boundary  of  P.  Note  that  since  wc  know  how  to  compute  a  triangulation  of  an 
N-gon  in  0(ATog  N )  time,  this  result  constitutes  a  significant  improvement. 

For  the  time  being,  we  will  assume  that  both  A  and  B  arc  vertices  of  P.  Wc  will  sec  later  on  how  we  can 
easily  dispense  with  this  requirement  If  A  and  B  arc  vertices  of  the  same  triangle  of  T.  it  is  clear  that 
IPf  A.B)=  AB.  so  wc  may  assume  that  this  is  not  the  case.  In  the  following,  wc  will  say  that  an  interior  edge  of 
7*  is  A  ^crossing  if  its  endpoints  «.»  are  such  that  A,u.B.v  appear  this  order  around  the  boundary  of  P.  Let  P' 
be  the  polygon  resulting  of  the  removal  from  T  of  ail  the  edges  that  arc  not  All-crossing  (fig.7).  Wc  first  prove 
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Proof)  It  suffices  to  show  that  IP tA,B)  can  only  intersect  AR-crossing  edges.  To  see  dial, 
suppose  that  it  intersects  an  interior  edge  e  which  is  not  AB-crossing.  Since  e  partitions  P  into  two 
polygons,  one  of  them  does  not  contain  B,  therefore  IFM.i)  crosses  e  at  least  twice  (once  in  each 
direction).  If  A  (resp.  B*)  is  the  first  (resp.  second)  intersection,  going  from  A  to  B,  replacing  the . 
pan  ofIPUB)  from  A *  to  B*  by  the  segment  a"B *  will  shorten  the  length  of  IPM.£),  which  leads 
to  a  contradiction.  □ 

Leona  4:  The  internal  path  between  A  and  B  intersects  every  interior  edge  of  P*  exactly  once, 
and  intersects  no  other  edge  in  T. 


I 


& 


!.> 


Proof)  The  proof  of  Lemma  3  shows  that  IP(A.B)  cannot  intersect  any  interior  edge  more  than 
once.  On  the  other  hand,  we  can  easily  prove  by  induction  that  since  every  interior  edge  of  P* 
partitions  this  polygon  into  two  parts,  neither  of  which  contains  both  A  and  B,  it  must  intersect 
IPM.#)  at  least  once.  Putting  this  result  together  with  Lemma  3  completes  the  proof.  □ 

It  is  easy  to  compute  P*  in  O (JV)  time.  To  do  so,  consider  every  interior  edge  of  T  in  turn,  and  if  it  is  not 
AB-crossing,  remove  it  from  faking  with  the  dangling  sub-polygon, just  created,  that  does  not  contain  A  at 
B.  Let  L  *  be  the  interior  edges  of  f*.  as  they  appear  from  A  to  £(fig.7).i.e,  in  the  order  in 

which  they  intersect  IPfi4.ll)  (Lemma  4).  Note  that  it  is  straightforward  to  obtain  L  in  0(A)  time,  once  P*  has 
been  computed.  From  now  on.  the  term  IFUofl.  with  x.y,  vertices  of  P*.  refers  to  the  internal  path  between  x 
and  y  with  respect  to  either  P  or  P*.  This  is  legitimate  since  the  two  paths  are  identical  as  a  simple 
generalization  of  Lemma  3  readily  shows. 

(FIGURE  8] 

Figure  8:  Computing  IP(A.B  iteratively.) 


H 


Lemma  5:  For  any  i;  lgisp,  there  exists  a  vertex  v  of  P*  such  that  IP(/f.^)s  IPfyf.v)  u  (/and 
IPM.tyz  IPM.v)  U  W,  where  U  and  IK  are  two  convex,  non-intersecting  polygonal  lines  mming 
their  convexity  against  each  other,  and  running  from  v  to  a.  and  respectively  (flg.8). 

Proof)  Let  C1  and  C,  be  two  oriented  curves  originating  at  the  same  point.  To  carry  the  analogy 
with  internal  paths,  we'may  forthcr  assume  that  neither  is  self-intersecting:  we  say  that  C,  and  C2 
have  a  proper  crossing  if.  as  we  follow  C,  from  its  suiting  point,  wc  encounter  a  point  where  C2 
intersects  Cr  and  actually  switches  from  one  side  to  the  other.  Fig.9.1  (but  not  fig.9.2)  shows  an 
example  of  a  proper  crossing. 

(FIGURE  9] 

Figure  9:  The  notion  of  proper  crossing. 
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We  next  prove  that  for  any  three  points  A.B.C  in  f.  the  two  paths  VHAB)  and  IF(AO  never 
have  any  proper  crossings.  Suppose  that  they  did:  let  a  be  the  first  point  (starting  at  A)  when 
IP {AB)  and  IP(  AO  cease  to  coincide,  and  let  b  denote  the  next  intersecting  point  Since 
and  IPlAO  take  distinct  paths  from  a  to  6.  we  may  re-route  either  one  to  die  other,  since  (hey 
must  have  exactly  the  same  length.  Iterating  on  this  process  win  eventually  cause  an  proper 
crossings  to  disappear,  which  proves  the  above  fact.  We  can  now  establish  Lemma  S  by  induction 
mi  L  The  initial  case  befog  trivial  we  may  directly  assume  that  the  lemma  is  true  for  all  indices 
from  1  to  i  Since  the  a^b-’s  are  triangulation-cdges.  we  necessarily  have  a.{  ,  or  say, 

wl°8*  Tbus>  considering  the  path  IPfAb^),  we  observe  that  since  it  does  not  have  any 
proper  crossings  with  either  IPMur,)  or  IPf/lb^ 

L  It  must  pass  through  their  common  point  r. 

1  Its  vertices  between  v  and  b^+1  are  vertices  of  U  and  W. 

From  L.  it  results  that  we  may  concentrate  on  the  path  IP(v,6i+1)  instead  of  IPM.h.  A  since  we 
obviously  have  IFLUl+.)  =  1P(  Av)  u  IP(v.bj+1).  Next,  wc  strengthen  proposition  1  by  proving 
that  die  vertices  of  IP(v.^1)  are  vertices  of  U  or  W.  but  never  of  both  at  the  same  time.  Indeed. 

suppose  wfog  that  starting  at  v.  the  vertices  of  IPC  are  ivu _ with  /j  through  t  lying  on  U 

and  fm+1  on  W.  It  follows  that  die  angle  1>r_/nH)  is  under  180  degrees,  therefore  there  is 
an  obvious  shortcut  for  IPfr,^  A  avoiding  (fig.  8)).  which  leads  to  a  contradiction.  Thus  there 
are  now  two  basic  cases  to  consider,  depending  on  whether  takes  its  vertices  in  (for  W. 

In  the  former  case,  »  will  be  relocated  ftirther  ahead  on  U.  whereas  it  win  stay  imrhmyri  fo  the 
Utter.  The  deuib  are  straightforward,  so  we  may  conskier  the  proofof  the  tanma  as  complete.  □ 

[FIGURE  10  J 

Rpn  10:  Minimality  propgrtitscflFfub^ilL 

We  are  now  ready  to  proceed  with  the  algorithm  for  computing  VP{ABA  The  method  in vohes  computing 
lP(^.u,)  and  IF(46jX  for  which  we  can  do  iteratively  by  using  the  reads  of  Lemma  S.  Tha 

procedure  befog  trivial  for  i*L  we  turn  to  the  general  step  directly.  As  already  mentioned,  we  have  either 
^  A  or  b|*  bi+  r  and  we  can  assume  wlog  that^*a,+r  Let  ul^Mm  (rasp.  wt — w^)  be  the  vertices  of  V 

(resp.  HO  from  v  to  a,  (resp.  bp. 

The  half-plane  delimited  by  afa on  (beside  where  4^+l  lies  is  partitioned  into  a+0+1  regions,  themselves 
delimited  by  the  lines  passing  through 

With  titis  order,  the  regions  appear  sorted  along  the  segment  a,/*  from  b(  to  aj,  so  that  wc  can  find  tha 
region  which  contains  b)+1  by  testing  each  of  them  in  turn  in  this  order,  until  we  are  successful  (ftg.ll).  This 
corresponds  to  unfolding  W  and  possibly  folding  over  U.  Ifb.+  1  lies  in  a  pencil  of  the  kind  (’Vi*V,W+ 
wc  must  simply  remove  from  W  and  reset  ft  to  k+ 1  and  to  ,  (fig.1 1.1).  If  bl+ ,  lies  in  tilt 

pencil  however,  wc  must  set  W  to  Vi*v  remove  (v.i/j _ «k  l)  from  V  and  finally  set  v  to 

(fig.1 1.2).  All  the  other  cases  are  similar  and  call  for  no  further  explanation.  Since  none  of  the  vertices 


removed  in  these  operations  will  ever  be  examined  again  (Lemma  5),  both  IP [AmJ  and  IPM.^  hence 
IPL4.A).  will  be  computed  in  0 (JV)  time. 

I  FIGURE  11] 

Figure  11:  Updating  Vtmd  W. 

We  generalize  this  result  by  allowing  both  A  and  B  to  Be  anywhere  inside  F,  and  not  only  on  the  boundary. 
Let  k  (resp.  5)  be  the  triangle  where  A  (reap.  B)  ties.  If  R»S,  the  problem  is  solved  since  YP(A,B)*AB. 
Otherwise.  we  can  compute  the  chain  of  triangles  /**  in  exactly  the  same  way  as  described  above. 

Next,  let  be  the  interior  edge  of  R  which  IP(iiJ)  cremes.  We  can  replaced  by  the  triangle  y^A  without 
altering  the  path  IF(.i4.R).  Applying  the  same  treatment  lo  5  win  make /<  and  £  become  vertices  of  P*.  which 
allows  us  to  call  on  the  procedure  described  eatttsr  to  compute  IP(X,£).  In  conclusion,  we  can  state  our  main 
resale 

Theorem  fc  Let  P be  a  staple  N-pm,  and  assume  that  any  triangulation  of  Fis  available.  For 
iuy  pair  of  points  A.Bm  F,  it  is  pomibie  recompute  IFf/LFL  the  unernal  path  between  A  and  B, 

3.2.  Thw  aM*intwmaH>attt«  praMuaa 

The  prablem  is  re  pnpracaas  the  polygon  Fro  that  a  batch  of  queries  of  the  Had: 

What  it  Bn  internal  pmtk  bttmm  A  emf  J? 

can  be  answered  optimally.  The  method  dcacriBcd  in  the  previous  section  grants  an  attractive  balance 
between  execution  and  preprocessing  time,  when  only  a  few  queries  have  to  be  handled  at  any  given  time.  It 
is  worst-case  optimal  but  not  optimal  in  the  strictest  sense  of  the  term.  since  all  the  vertices  of/*  must  always 
be  examined  fcr  every  query.  As  a  result,  the  precomputation  of  all  possible  internal  paths  between  vertices 
entails  a  prohibitive  OUV8)  cost  The  goal  which  we  set  forth  here  is  re  pwproccsi/*  so  that  the  computation  of 
IF {A.B)  hr  any  pair  of  vertices  (A.B)  requires  only  time  proportional  to  the  size  of  the  output.  Le*  the 
number  of  vertices  in  IF(i4,£). 

To  achieve  this  goal  we  use  the  concept  of  visibility  introduced  earlier.  Let  V(A)  be  the  visibility  polygon 
of  A.  If  IFM.£)m  AB.B  is  a  vertex  of  V(^),  otherwise  V(A)  has  a  ray-edge  (i.e.  an  edge  vwsuch  that  v  lies  on 
Aw).  with  the  property  that  vw  separates  A  from  B  by  intersecting  IPM.R).  More  precisely,  vw  is  the  unique 
edge  of  V(/|)  such  that  either  A,y.B,w  or  A.w.B,v  occur  in  clockwise  order  (fig.12).  Since  V(/t)  is  star-shaped, 
and  vw  is  a  ray-edge  which  is  traversed  by  IP(  A.B).  v  must  be  the  first  vertex  of  IP(  A.B)  after  A.  Indeed,  there 
would  be  a  shortcut  if  1P(  A.B)  cut  vw  at  any  other  point  Consequently,  we  have  the  relation: 

IP(/f,£) »  /tvu  IF(v,£) 
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This  motivates  the  introduction  of  the  function  F {A,B)~  B,  if  IPM,5)=  AB,  and  F{A,B)=  v  otherwise. 
Theorem  2  shows  that  if  a  triangulation  of  P  is  available,  the  visibility  polygon  V(/f)  of  each  vertex  A  of  Fean 
be  obtained  in  0(N)  time.  The  knowledge  of  V(/f)  permits  us  to  set  up  the  army 

D^O*{FUVj);i=UJV} 
in  0(N)  tune,  with  0(N)  storage,  from  which  we  conclude: 

Theorem  7:  Let  P  be  a  simple  polygon  with  N  vertices.  It  is  possible  to  preprocess  P  in  OfM2) 
time,  using  CXJV2)  space,  so  that  for  any  pair  of  vertices  A,B,  the  path  IP (A,B)  can  be  computed 
optimally,  Le.  in  time  proportional  to  the  size  of  the  output 

Prooft  Compute  the  N  arrays  D(Vj)^D(v  J,  forming  an  NxN  matrix  {Ffv.v.}.  so  that  IP(A,B) 
can  be  computed  by  retrieving  r (A,3)  in  constant  time,  and  computing 
lP(F(/t,F),F)  recursively.  □ 

[FIGURE  12] 

Figure  12:  The  alF  internal- paths  problem. 


3.3.  The  intomaMongth  problem 

Imagine  that  an  island  with  only  inland  communications  is  to  be  serviced  by  some  utility  (water  tank, 
power  station,  file  house,  police  station,  hospital.  esc~).  An  interesting  piece  of  information  which  may  be 
needed  is  an  upper  bound  on  the  internal  path  length  between  any  pair  of  points. 

Let  A*.b‘  be  the  two  vertices  of  P  which  form  the  longest  path  IP {A*.B*).  We  call  |IP(/I*.F*)|  the  intend 
length  of  P.  It  is  easy  to  determine  A  *  and  B'  by  trying  out  all  possible  pairs  of  vertices  and  using  the  matrix  F 
of  the  previous  section,  given  that  tire  longest  path  can  always  be  assumed  to  be  found  between  two  vertices  of 
foe  polygon.  This  leads  toanCXAf5)  running  time,  which  we  can  cut  down  to  O (fV2)  by  proceeding  as  follows: 

Let  D (A.BJ  =  |IP(/f.F)J.  We  will  compute  D (A.B)  iteratively  by  summing  up  partial  distances  obtained 
from  F.  In  order  to  avoid  duplicating  computations,  as  soon  as  TXA.B)  is  available,  wc  backtrack  along  foe 
path  just  followed  in  F  to  record  the  partial  results.  This  ensures  that,  on  average,  one  value  D (A,B)  will  be 
computed  at  every  other  step,  which  leads  to  an  O (N2)  algorithm. 
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INTDIST 


-  Initially,  each  D {A,B)  is  set  to  -1  for  A*B,  and  to  0  for  A=B. 

for  all  i  (l^isiV) 
for  all  jfl^jsAO 

begin 

<2— <v,} 

while  D(jf.v.)  =  >1 


at<-F(x.vp 

end 

if  Q  has  more  than  one  element 


Let  Q  = 

I— IX^vJ 
for  k  =  p-I^.,1 


Z.«~L+|^+1| 

D(J 


D {A*.B*)  =  Max  ( D(  *,  vp  |  all  pain  of  vertices  vp  Vj ) 


Since  wc  can  compute  a  triangulation  of  Fin  OfA/log  AO  time,  we  may  conclude: 

Theorem  8:  It  is  possible  to  determine  the  internal  length  of  a  simple  Af-gon  as  well  as  the 
corresponding  internal  path  in  O (A/*)  time. 


4.  Conclusions 

This  paper  has  shown  on  the  following  examples  how  to  use  an  arbitrary  triangulation  advantageously: 

1.  Computing  the  visibility  polygon  at  any  point  inside  an  N-gon  in  0(N)  time. 

2.  Computing  the  internal  path  between  any  pair  of  points  in  an  N-gon  in  0^^ 0  time. 

3.  Allowing  0(N2)  preprocessing,  being  able  to  compute  any  internal  path  with  optimal 
performance. 

4.  Computing  the  internal  distance  of  an  N-gon  and  the  associated  internal  path  in  0(Af2)  time. 
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All  of  these  algorithms  achieve  significant  improvements  over  previously  known  methods,  since  a 
triangulation  of  an  N-gon  can  be  computed  in  0(Mog  N )  time.  The  improvements  are  to  be  measured  either 
in  terms  of  better  performance  (Problems  23.4)  or  in  terms  of  added  simplicity  (Problem  1).  We  should  also 
observe  that  it  is  yet  unknown  whether  the  triangulation  algorithms  available  in  the  literature  are  optimal. 
Since,  on  the  other  hand,  half  of  the  algorithms  which  we  have  described  in  this  paper  are  linear  after 
triangulation,  overall  speed-ups  would  automatically  result  from  the  discovery  of  hater  triangulating 
procedures. 

This  work  was  meant  as  a  case-study  and,  of  course,  the  list  of  possible  improvements  brought  about  by  the 
use  of  a  triangulation  is  not  dosed.  Further  research  should  attempt  to  enlarge  the  list  given  here,  and  cany 
the  same  approach  with  other  preprocessing  structures,  whether  geometrical  or  not 
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